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Abstract. We study the large-scale angular correlation signatures of the Cosmic Microwave Background (CMB) 
temperature fluctuations from WMAP data in several spherical cap regions of the celestial sphere, outside the 
KpO or Kp2 cut-sky masks. We applied a recently proposed method to CMB temperature maps, which permits an 
accurate analysis of their angular correlations in the celestial sphere through the use of normalized histograms of 
the number of pairs of such objects with a given angular separation versus their angular separation. The method 
allows for a better comparison of the results from observational data with the expected CMB angular correlations 
of a statistically isotropic Universe, computed from Monte Carlo maps according to the WMAP best-fit ACDM 
model. We found that the, already known, anomalous lack of large-scale power in full-sky CMB maps are mainly 
due to missing angular correlations of quadrupole-like signature. This result is robust with respect to frequency 
CMB maps and cut-sky masks. Moreover, we also confirm previous results regarding the unevenly distribution 
in the sky of the large-scale power of WMAP data. In a bin-to-bin correlations analyses, measured by the full 
covariance matrix y 2 statistic, we found that the angular correlations signatures in opposite Galactic hemispheres 
are anamalous at the 98%-99% confidence level. 
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1. Introduction 

In recent years, after the COBE mission, experiments 
of observational cosmology evolved both in the accuracy 
and in the angular resolution of the Cosmic Microwave 
Background (CMB) measurements (see, e.g. Bersanelli et 
al. 120021 and references therein). However, such experi- 
ments were designed to map small patches of the celestial 
sphere and therefore the study of its properties concerning 
the large-scale angular correlations was not possible. 

This situation changed with the Wilkinson Microwave 
Anisotropy Probe (WMAP), which produced full-sky 
CMB maps in five frequency bands (termed K, Ka, Q, 
V, and W maps), although containing different amounts 
of contaminations from our galaxy (Bennett et al.[2003a). 
This superb quality data set motivated a number of de- 
tailed studies concerning the CMB properties. As a re- 
sult, several recent works have reported claims regarding 
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some anomalies found in these data. These anomalies in- 
clude the low-order multipole values (Bennett et al. 2003a 
Efstathiou |2CjOH Tegmark et al. l2~0"0"3l Gaztahaga et 
al. I2UU51 Slozar et al. I2UU1 O'Dwyer et al. EUUljl . 
the alignment of some low-order multipoles (Tegmark 
et al. 120031 de Oliveira-Costa et al. 2004, Eriksen et 
al. l2004bl Bielewicz et aL lStMl Land and Magueiio lMTSjl . 
and an unexpected asymmetric distribution on the sky 
of the large-scale power of CMB data (Eriksen et 
al. l2^04all2UUoal Hansen et al. I^004al I2004E1 l^004cl Jaffe 
et al. l2Utj51 Schwarz et al. lSUUl Copi et al. lSUUfjjl . 



In this work we focus on the study of the large-scale an- 
gular correlations (hereafter termed Angular Correlations 
Signatures ACS) of the CMB temperature maps from 
WMAP data. Since the vast majority of the cosmologi- 
cal information is contained in the two-point temperature 
correlation function, it seems natural to study the ACS 
by looking at such function (actually, an equivalent of it). 
Our analyses of the ACS in CMB maps, for data outside 
the region defined by the KpO or Kp2 WMAP masks, al- 
low to perform a close inspection of two interesting issues. 
First, we investigate -through a full-sky map analysis- 
the large-scale power in WMAP data, and its connection 
with the low quadrupole moment issue. Second, we study 
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-through a partial-sky coverage analysis- the possible un- 
even distribution on the celestial sphere of the CMB large- 
scale power. The significance of our results are evaluated 
by comparing these results against those obtained from 
Monte Carlo CMB maps produced with the WMAP best- 
fitting ACDM model properties (Hinshaw et al. 12003b). 
Finally, using a \ 2 statistics we assess the confidence level 
of the North/South asymmetry in WMAP data compared 
to 1 000 Monte Carlo realizations of CMB skies. 

2. Analyses method 

Recently, Bernui & Villela (2006) introduced a new 
method to study the large-scale ACS in the distribution 
of cosmic objects in the sky. This method, called the Pair 
Angular Separation Histogram (PASH) method, consists 
in first calculating the angular distances between all pairs 
of cosmic objects, listed in a catalog, and then construct- 
ing the normalized histogram of the number of pairs with 
a given angular separation versus their angular separation. 
A catalog with a large number of objects can be divided 
in a number (say K) of comparable sub-catalogs (i.e. en- 
sembles of iso-number of objects sharing analogous phys- 
ical properties). After that, one computes a PASH with 
each one of these sub-catalogs and average them to ob- 
tain the Mean-PASH (MPASH). The difference between 
the MPASH, caculated using a observational catalog, and 
the Expected- PASH (EPASH) , obtained assuming the sta- 
tistical isotropy hypothesis (hereafter the MPASH-minus- 
EPASH function), shows the ACS of the cosmic objects in 
such a catalog. 

In the CMB temperature maps the celestial sphere is 
partitioned in a set of equal-area pixels, where to each 
pixel is assigned a weighted CMB temperature. To obtain 
the ACS of a given map, one divides the set of pixels in 
K = 2 sub-maps, one for the negative CMB temperature 
fluctuations and the other for the positive ones, and pro- 
ceeds as before averaging (K —) 2 PASHs. If, for compu- 
tational problems, the number of pixels in these sub-maps 
is too large, one can divide them in, say A_ and K + , 
sub-maps (with iso-number of pixels of analogous tem- 
peratures), and proceeds to compute the MPASH aver- 
aging the K = K- + K + PASHs, and finally, one plots 
the MPASH-minus-EPASH function. It can be shown (see 
Bernui IMiS) that the MPASH-minus-EPASH function is 
independent of the number of sub-maps A"_ > I and 
K + > 1 in which the original CMB map is divided. 

The PASH method is also applicable to incomplete 
sky maps, including disconnected regions such as those 
resulting from the application of the WMAP masks. This 
method is similar in philosophy to the two-point temper- 
ature correlation function, except that the PASH method 
has zero mean, because the MPASH and EPASH are each 
one normalized histograms. This simple fact allows us to 
achieve a deeper study of the ACS in CMB temperature 
maps, like: the intensity of all their relative mimima and 
maxima and their corresponding angular scales, their in- 
tersections with the horizontal axis (i.e. the zeros of the 



function), etc; it also leads to calculate the variance a 2 
of the MPASH-minus-EPASH function ft (which has zero 
mean, i.e., ^=1" ft = °) 

2—1 

3. CMB data analyses 

In this section we present the ACS of the WMAP data. 
For these analyses, we used six CMB maps: the individual 
Q-, V-, and W-band data produced from the (Qi,Qj2), 
(Vi,V 2 ), and (W1W2W3, W4) differential assemblies 
(DA), respectively, as provided on LAMBDA 1 (Bennett 
et al. 12003 bjl . the Coadded WMAP map (which com- 
bines the eight DA in the Q-,V-,and W-bands using the 
inverse- variance noise weights as described by Hinshaw et 
al. I|2003bjl ). the cleaned CMB map (Tegmark et al. l2"0"0H 
hereafter the TOH map), and the Lagrange Internal 
Linear Combination (Eriksen et al. I2004bl hereafter the 
LILC map). It is known that the last two maps are not 
fully reliable for quantitative analysis because still have 
problems with residual foregrounds (Eriksen et al. 1200 5b): 
here we use them just as a supporting evidence. After 
applying the KpO or Kp2 WMAP masks in all these six 
maps, we remove their residual monopole and dipole com- 
ponents, and correct them for the dynamic quadrupolc. 

The examination of the WMAP data performed here 
comprehends two types of sky analyses. We want to remark 
that in both cases we consider data outside the region de- 
fined by the KpO or the Kp2 mask, to avoid from the 
beginning residual foreground contaminations (Bennett 
et al. I2003bl Eriksen et al. I2004bl I2005bl Bielewicz et 
al. I2004|) . The first type is the full-sky analysis, and the 
second one is the partial-sky analysis carry out through 
antipodal spherical caps of 45°, 60°, and 90° of aber- 
ture (this means that the maximum angular separation 
between pixels in each case is twice these angle values). 

We use CMB maps with HEALPix (Gorski et al. 
resolution N si d c = 128, which amounts to 196 608 pixels. 
In the following figures we consider the bin width as 0.45°, 
and the number of bins as JVbi ns = 400; moreover, the 
shaded areas correspond to 1-sigma confidence regions. 

3.1. Full-sky analyses (outside Kp2 or KpO masks) 

In this subsection we shall discuss the ACS of the WMAP 
maps, analysing the CMB data outside the regions defined 
by the Kp2 and KpO masks. The objective of such anal- 
ysis is to investigate previous claims regarding the lack 
of power in the large-scale angular correlations in full-sky 
data (see Bennett et al. B003a Efstathiou 120041 Slosar et 
al. I2004J1 and its relationship to the low-order multipole 
values. 

1 http:\\lambda.gsfc.nasa.gov\product\map\current 
\IMaps_ cleaned, cfm 
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First we show the ACS of the WMAP data using the 
Kp2 mask. In fact, the suitability of using each WMAP 
mask is well known in the literature, and many works 
claim that the less severe cut-sky Kp2 could be enough 
for cosmological purposes (large angular scales). We have 
computed the ACS for both masks, although we only show 
the ACS for the Kp2 mask case, and for all maps under 
investigation to assert the robustness of our results under 
different sky cuts. 

The results are reported in Figure ^ In the top panel 
of Fig.^we show the ACS from the six CMB maps here in- 
vestigated, considering data outside the region determined 
by the Kp2 mask. For comparison we also plot, as a dashed 
line, the average of 1 000 MPASH-minus-EPASH functions 
(termed the expected function) obtained from the same 
number of Monte Carlo CMB maps. These realizations 
were produced using an input power spectrum generated 
by the CMBFAST code (Seljak & Zaldarriaea I1996[l con- 
sidering the WMAP best-fitting ACDM model properties. 
The maps were produced using the map-making code SYN- 
fast from HEALPix. The disagreement observed between 
the observational data from WMAP and the expected an- 
gular correlation function is better appreciated in the bot- 
tom panel of Fig. ^ where we plotted the difference of 
such expected function minus each one of the MPASH- 
minus-EPASH functions from the six CMB maps. What 
is revealed in these difference plots is that the ACS miss- 
ing in CMB data at large-angular scales correspond to 
a quadrupole-like angular correlation signature (see, e.g., 
Bernui 2005 ). To assess this, we plot the angular correla- 
tion function of a quadrupole (corrected for the Kp2 mask) 
with moment C*2 — 670/iK 2 . Since these curves (bottom 
panel Fig.Q represent the missing ACS in WMAP data to 
fit the expected function, where C^ CDM = 870/iK 2 , one 
finds Cf MAP ~ 190 ± 30^K 2 , that are in good agree- 
ment with Bielewicz et al. J2004) who obtained C*2 = 
165 ± 34, 216 ± 42, 229 ± 44, 191 ± 38 /jK 2 , for the Q, V, W, 
and Coadded maps, respectively). This result, which ap- 
pears to be robust with respect to the frequency CMB map 
and to cut-sky mask (we also performed analyses with the 
KpO mask obtaining similar results), proves that the lack 
of angular power previously found in the two-point tem- 
perature correlation function of WMAP data is mainly 
due to the low quadrupole moment (Efstathiou 2004 
Gaztanaga et al. 120031 Slosar et al. I2004j> . 

Furthermore, to confirm this result we remove the 
quadrupole component, after applying the Kp2 mask, and 
compute the ACS in these maps. The result is shown 
in Figure [21 where we also plot for comparison the ex- 
pected function for the case of zero quadrupole compo- 
ment. Again this result is robust with respect to the fre- 
quency CMB map analyzed. 

3.2. Partial-sky analyses 

The basic idea is to investigate different regions of the ce- 
lestial sphere, namely spherical caps of various sizes, whose 
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Fig. 1. Top panel: ACS in the six CMB maps: Q (violet 
line), V (blue line), W (light-blue line), Coadded (green 
line), TOH (brown line), and LILC (red line), after apply- 
ing the Kp2 mask. For comparison we also plot the av- 
erage of 1 000 MPASH-minus-EPASH functions obtained 
from the same number of Monte Carlo CMB maps (dashed 
line, also called the expected function). One unit in the 
vertical axis corresponds to 5 460/iK 2 . Bottom panel: The 
differences of the expected function minus the six MPASH- 
minus-EPASH functions from the top panel, where we ob- 
serve that these ACS correspond mainly to a quadrupole- 
like signature. The dashed line here corresponds to the 
MPASH-minus-EPASH of a quadrupole (corrected for the 
Kp2 mask) with C2 — 670/iK 2 . Notice that these curves 
represent the missing ACS in WMAP data to fit the ex- 
pected function, which have Cf CDM = 870/iK 2 . 

vertices point in the direction of the North Galactic Pole 
(NGP) and the South Galactic Pole (SGP), with respect 
to the Galactic coordinate system. Comparison of the ACS 
in antipodal caps of 45°, 60°, and 90° of aberture allows 
a quantitative asscsscmcnt of the claimed asymmetry in 
the distribution of the CMB power between the northern 
and southern hemispheres (Eriksen et al. l2004Ej) . 

The resulting ACS are shown in Figures|3and0J where 
the expected functions result from data on the same patch 
of the sky from Monte Carlo CMB maps. In Figure [31 
we present four plots: panels (a) and (b) reveal the ACS 
in the 45° NGP-cap and SGP-cap, respectively; while in 
panels (c) and (d) we show the ACS in the 60° NGP- 
cap and SGP-cap, respectively. In Figure 21 we also show 
four plots: panels (a) and (b) reveal the ACS in the 90°- 
cap, outside the Kp2 mask, in the Northern and Southern 
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Fig. 2. The ACS for the six CMB maps (the nomenclature 
of the curves follows the same color's pattern as before) 
with the quadrupole component removed. The dashed line 
is the expected function for the case of zero quadrupole 
moment. 

Galactic hemispheres (NGH, SGH), respectively; panels 
(c) and (d) show the ACS in the 90°-cap, outside the Kp2 
mask, in the NGH and SGH, respectively, but this time 
we removed the quadrupole component of the CMB maps. 

Some interesting features in these plots are evident. 
The first noticeable result is that, in each one of these fig- 
ures the ACS corresponding to the three frequency maps 
Q, V, and W are practically indistinguishable. It is well 
known (see Bennett et al. 2003b|) that foregrounds have 
a frequency-band dependence, thus when incorrectly sub- 
tracted it is expected to show up at a different level in 
the frequency maps Q, V, and W. To investigate the ro- 
bustness of our results with respect to a different sky cut, 
we computed the ACS for the NGH and SGH (not shown 
here) for data outside the more severe galactic cut given 
by the WMAP KpO mask (it cuts ~ 25% of the total 
data, while the Kp2 mask cuts ~ 15%), where we found 
similar results. Thus, one concludes that the ACS found 
in these maps are remarkably stable with respect to both 
frequency bands and sky cuts, and seem unlikely to be 
compromised by residual foregrounds. 

The second interesting fact is the net asymmetry be- 
tween the 90°-caps in NGH and SGH: the Northern hemi- 
sphere is almost structureless (the variance's square-root 
for this ACS is (cr$' V ' W ) = 0.0054 while ^Expected = 
0.012 ± 0.006), instead the ACS in the Southern hemi- 
sphere are relatively intense ((<7g ^'^) = 0.016 while 
^Expected = o.015±0.008). In this form, the South/North 
ratio (?s/n = vs/vn = 2.97 quantifies the evident 
SGH/NGH asymmetry observed in the plots of figure 0] 
The statistical significance of this result are evaluated in 
the next section, by means of the Monte Carlo maps, in 
terms of the standard covariance matrix \ 2 statistic. 

Finally, in order to check the possibility that the low 
value of the quadrupole moment could be responsible for 
this peculiar ACS behavior, we removed the quadrupole 
component, after applying the Kp2 mask, in all the 
six maps, and then compute again the MPASH-minus- 
EPASH for data in the 90°-caps in the NGH and SGH. 



Our results can be seen in the Figures andQJl, where 

this time (cr$' V ' W ) = 0.0065 in the NGH, and (ctJ' V,W ) 
= 0.019 in the SGH. The large value now obtained for 
the ratio <Js/n — 2.96, confirms that the North/South 
asymmetry is not related to the low quadrupole value. 

A natural question concerns the frequency that such 
a ratio is os/n — 3 {<Js/n — 3.2 if data from Coadded 
map is included in the average) is present in the ensem- 
ble of Monte Carlo maps. Analysing the MPASH-minus- 
EPASH functions for similar NGH and SGH 90°-caps in 
1 000 Monte Carlo maps, we found that the corresponding 
values of a s /N £ [0.2, 5.4], with mean aff$ n = 1.1 ± 0.7. 
A close inspection of the values in such interval results in 
that less than 2% of them satisfies os/n > 3. 

4. Statistical analyses 

In this section we present the statistical analysis of the 
ACS obtained from WMAP data against analogous re- 
sults derived from the 1 000 Monte Carlo CMB maps. The 
degree of agreement between the simulations and the ob- 
servations are quantified in terms of a standard covariance 
matrix \ 2 statistic, including all bin-to-bin correlations, 

iV binE , 

where is the mean of the MPASH-minus-EPASH 
functions as determined by the Monte Carlo maps, and 
Mij = (fifj) — (fi)(fj) is the covariance matrix. The re- 
sults from these computations, divided in three groups for 
the 45°-, 60°-, and 90°-caps, respectively, are reported in 
Table 1. The first column indicates the CMB map investi- 
gated, the second and third ones indicate the frequencies 
of Monte Carlo simulations with a lower \ 2 value than 
the WMAP data for the NGH and SGH, respectively. The 
last column indicates the frequency of simulations with a 
smaller XgGp/xJjQp- These results determine a parame- 
ter, the ratio of the \ 2 values for the NGH and SGH, that 
is used to quantify the degree of asymmetry between both 
hemispheres (Eriksen et al. l2004al l2~005al. 

The results reported in Table 1 are as follows. In all the 
cases, the NGH and SGH show, indepedently, acceptable 
values of probability (frequency of simulated maps, given 
the model, where the MPASH-minus-EPASH function of 
such maps has a lower \ 2 value than the WMAP data). 
However, data in the third column says that the antipodal 
90°-caps are marginally consistent internally (when con- 
sidered simultaneously one as opposite to the other): for 
the Q, V, W, and Coadded maps the mean of the x%/n 
values is ^ 0.98, or equivalently X%/s & 0.02, which means 
that 5s 2% of the simulations have a smaller ratio of \ 2 - 
It is worth to notice that these results appearing in the 
90°-caps are shown to be robust with respect to frequency 
(the same result for Q, V, W, and Coadded CMB maps). 
Therefore, we conclude that there is a significant discrep- 
ancy, at the ^ 98% CL, between WMAP data and the 
ACS expected in a statistically isotropic Universe. 
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Table 1. x 2 -test for the ACS in different caps. The first 
column indicates the CMB map investigated, the second 
and third ones indicate the frequency of Monte Carlo sim- 
ulations with a lower % 2 value than the WMAP data. The 
last column indicates the frequency of simulations with a 
smaller x| GH /XnGH' 



CMB map 




X 2 -NGH 


X 2 -SGH 


2 

Xs/N 




45° 













0.914 


0.891 


0.499 


V 




0.854 


0.902 


0.624 


w 




0.654 


0.782 


0.616 


COADDED 




0.969 


0.967 


0.548 


TOH 




0.865 


0.759 


0.419 


LILC 




0.817 


0.882 


0.626 




60° 













0.715 


0.971 


0.949 


V 




0.637 


0.982 


0.956 


w 




0.361 


0.959 


0.930 


COADDED 




0.550 


0.909 


0.822 


TOH 




0.786 


0.926 


0.712 


LILC 




0.154 


0.884 


0.891 




90° 








Q 




0.276 


0.959 


0.985 


V 




0.243 


0.810 


0.965 


w 




0.278 


0.968 


0.976 


COADDED 




0.143 


0.929 


0.990 


TOH 




0.164 


0.724 


0.846 


LILC 




0.176 


0.770 


0.853 



5. Conclusions 

Here we have shown through the full-sky analyses of the 
ACS of the WMAP data (outside the KpO or Kp2 sky 
cuts) that the 'anomalous' lack of power at large angular 
scales is mainly explained by the low quadrupole-moment 
value (see Efstathiou 200 1 Gaztahaga et al. 120031 and 
Slosar et al. 120041 for more detailed discussions; see also 
Bielcwicz et al. 2004 for complete calculations, in partic- 
ular in similar band maps as here). We confirmed this 
result by analysing also the WMAP data after removing 
the quadrupole component, obtaining a better (but not ex- 
actly) agreement with what is expected in a ACDM model 
with zero quadrupole. 

We also performed partial-sky coverage analyses and 
showed that there is significative evidence, at the 98% — 
99% CL, of an anomalous distribution of the CMB power 
in opposite Galactic hemispheres. In fact, in less than 2% 
of the Monte Carlo CMB maps analyzed, data in oppo- 
site hemispheres present analogous bin-to-bin correlations 
as those shown by the 90°-caps with WMAP data. More 
interestingly, the excess of correlations in the comparison 
of opposite antipodal caps is manifestly observed in the 
South/North 90°-caps, but it is not seen neither in the 
examination of the antipodal 45°- or 60°-caps. 

The fact that the South/North asymmetry persists at a 
significant level in several frequency bands (the Q, V, and 
W CMB maps) and with different sky cuts (the KpO and 
Kp2 masks) excludes a simple explanation based on sys- 



tematic effects and foregrounds (see Hinshaw et al. 2003a 
Eriksen et al. I2004al l2ljfJ4bl Hansen et al. l2004bjl . 
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Fig. 3. The ACS in the 45°-caps in the NGH (a) and in 
the SGH (b), and the ACS in the 60°-caps in the NGH (c) 
and in the SGH (d), from the six CMB maps (Q, V, W, 
Coadded, TOH, and LILC), after applying the Kp2 mask. 
In all plots the dashed line is the average of MPASH- 
minus-EPASH functions obtained from the corresponding 
sky regions in 1 000 Monte Carlo maps, respectively. The 
nomenclature of the curves follows the same color's pat- 
tern as in Figure n 
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Fig. 4. The ACS in the 90°-caps in the NGH (a) and in 
the SGH (b), from the six CMB maps (Q, V, W, Coadded, 
TOH, and LILC), after applying the Kp2 mask. Removing 
the quadrupole component, after applying the Kp2 mask, 
in the six CMB maps we obtain the ACS in the 90°-caps 
in the NGH (c), and in the SGH (d). For comparison, 
the dashed line and the shaded region in panel (a) (panel 
(b)) is the same as in panel (c) (panel (d)). As before, the 
nomenclature of all the curves follows the same color's 
pattern as in Figure Q 



